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Self-organized  quantum  dots  pattern  depends  strongly  on  the  elastic  strain  energy  of  the  substrate. 

It  is  well-known  experimentally  that  for  the  elastic  substrate  with  a  high  degree  of  anisotropy,  the 
epitaxially  grown  island  patterns  are  different  for  different  growth  orientations.  In  this  paper,  by 
incorporating  the  anisotropic  strain  energy  held  into  a  kinetic  Monte  Carlo  algorithm  for  adatom 
diffusion,  we  show  that  the  self-organized  island  pattern  on  the  surface  of  an  anisotropic  substrate 
is  closely  correlated  to  the  elastic  energy  distribution  on  the  surface.  The  anisotropic  substrates 
studied  are  GaAs  with  different  growth  orientations  (001),  (111),  and  (113).  An  isotropic  substrate 
Iso  (001),  reduced  from  GaAs,  is  also  investigated  for  the  purpose  of  comparison.  The  island 
patterns  on  these  substrates  with  and  without  elastic  strain  energy  are  presented.  Besides  the  effect 
of  substrate  anisotropy,  different  growth  parameters,  including  temperature,  coverage,  and 
interruption  time,  are  further  investigated  to  identify  the  optimal  growth  values.  It  is  observed  that 
the  strain  energy  held  in  the  substrate  is  the  key  factor  that  controls  the  island  pattern,  and  that  the 
latter  is  closely  correlated  to  the  substrate  orientation  (anisotropy).  Our  simulated  patterns  are  also 
in  qualitative  agreement  with  recent  experimental  growth  results.  ©  2006  American  Institute  of 
Physics.  [DOI:  10.1063/1.2213153] 


I.  INTRODUCTION 

The  fabrication  of  semiconductor  nanostructures  has  at¬ 
tracted  considerable  interest  in  recent  years  due  to  their 
wide-range  applications  in  optoelectronics  and  semiconduc¬ 
tor  devices.  Among  many  other  techniques,  the  direct  growth 
of  nanostructures  has  evolved  as  a  promising  approach  to 
achieve  devices.  Usually,  the  thin  him  layer  has  different 
material  property  as  compared  to  the  substrate  material,  and 
thus  such  structures  are  generally  known  as  heterostructures. 
Due  to  the  difference  in  material  property  and  hence  the 
difference  in  the  lattice  parameter,  lattice-misht  induced 
strains  occur  and  thus  the  thin  layer  deposited  on  the  sub¬ 
strate  will  break  up  into  coherently  strained  islands  to  release 
some  of  the  strain  energy.  During  epitaxial  growth,  a  strained 
him  can  relax  its  strain  by  one  of  the  following  three  mecha¬ 
nisms:  (1)  formation  of  dislocation,  (2)  surface  roughening, 
and  (3)  cracking.  The  second  mechanism,  also  called  the 
Stranski-Krastanov  growth  mode,  has  attracted  signihcant  at¬ 
tention  due  to  its  potential  as  a  technique  to  manufacture 
self-organized  quantum-dot  (QD)  and  quantum-wire 
structures.4  The  driving  mechanism  for  this  growth  mode 
transition  is  the  very  efficient  strain  energy  relaxation. 

In  self-organized  QD  island  research,  both  theoretical5,6 
and  experimental  studies  have  shown  that  the  interaction 
among  islands  via  their  elastic  strain  helds 13,14  may  lead  to  a 
gradual  improvement  in  size  homogeneity,  as  well  as  to  a 
more  uniform  lateral  island  spacing.  Holy  et  al.  have  further 
shown  that  the  vertical  and  lateral  correlations  in  self- 
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organized  QD  superlattices  can  be  explained  by  taking  into 
account  the  elastic  anisotropy  of  the  materials.10  In  these 
superlattices,  dot  correlations  are  induced  by  the  interaction 
among  dots  via  their  elastic  strain  helds.  Above  the  buried 
islands,  the  elastic  energy  distribution  on  the  surface  exhibits 
pronounced  minima  and  maxima  in  the  lateral  directions. 
Many  continuum  mechanics  approaches,  including  the  finite 
element  analysis  and  Fourier  transformation  method,  were 
proposed  to  study  the  effect  of  elastic  anisotropy  on  the 
lattice-misht  induced  strain  energy  held.15-17  While  the  con¬ 
tinuum  mechanics  approaches  reveal  various  interesting  fea¬ 
tures  showing  the  direct  role  of  elastic  anisotropy  on  the 
self-organized  growth  pattern,17,18  they  fail  to  connect  the 
growth  ordering  and  pattern  to  the  growth  parameters,  a  se¬ 
rious  drawback  in  the  continuum  mechanics  analysis.  Fur¬ 
thermore,  the  identification  of  a  set  of  optimal  growth  param¬ 
eters  would  be  most  useful  to  experimentalists. 

The  kinetic  Monte  Carlo  (KMC)  has  been  proposed  re¬ 
cently  to  study  QD  island  self-organization  by  many 
researchers 19-24  where  the  growth  parameters  are  directly  in¬ 
volved  in  the  simulation.25  It  has  been  shown  that  the  main 
growth  parameters  affecting  QD  island  distribution  patterns 
are  the  temperature  T,  flux  rate  F  to  the  surface  during  depo¬ 
sition,  surface  coverage  c,  and  growth  interruption  time 
^  21,25  gase(j  on  a  proposed  coupled  KMC,  the  authors  simu¬ 
lated  the  island  ordering  and  narrow  size  distribution  in  two 
dimensions  and  further  identified  a  set  of  optimal  growth 
parameters  for  the  assumed  model  system.25  However,  how 
the  substrate  anisotropy  (or  orientation)  affects  the  island 
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FIG.  1 .  Contours  of  the  normalized  elastic  strain  energy 
distribution  on  the  substrate  surface  of  isotropic 
Iso(001)  (a),  and  anisotropic  GaAs(001)  (b),  GaAs(lll) 
(c),  and  GaAs  (113)  (d). 


pattern  has  not  been  studied  so  far  within  the  coupled  KMC 
algorithm  where  the  long-range  strain  energy  field  is  in¬ 
cluded. 

Therefore,  in  this  paper,  we  apply  the  coupled  KMC 
self-organization  model  to  study  the  effect  of  elastic  aniso¬ 
tropy  on  the  QD  island  lateral  ordering  and  pattern  [a  two- 
dimensional  (2D)  model  only].  In  our  simulation,  the  atom 
hopping  rate  is  governed  by  the  Arrenius  law  which  is  en¬ 
hanced  by  the  long-range  strain  energy  field  due  to  the 
lattice-misfit  strain.  The  induced  elastic  strain  field  is  ob¬ 
tained  using  the  Green’s  function  solution  for  the  anisotropic 
semiconductor  substrate  developed  recently.26,27  Besides  the 
isotropic  substrate  model  Iso  (001),  three  anisotropic  sub¬ 
strates  derived  from  GaAs  were  included,  i.e.,  GaAs  (001), 
GaAs  (111)  and  GaAs  (113).  Different  growth  parameters, 
namely,  temperature  T ,  coverage  c,  and  interruption  time  th 
are  also  discussed.  While  results  with  different  growth  pa¬ 
rameters  demonstrate  that  temperature,  coverage,  and  inter¬ 
ruption  time  affect  the  average  island  size  and  number  of 
atoms  within  each  island,  our  simulations  further  clearly  in¬ 
dicate  the  close  correlation  between  the  substrate  orientation 
(elastic  anisotropy)  and  the  growth  pattern,  which  also  quali¬ 
tatively  agree  with  recent  experimental  results  obtained  via 
scanning  tunneling  microscopy  (STM).28,29 

II.  SELF-ORGANIZED  ORDERING  AND  PATTERN 

In  our  coupled  KMC  simulation,  the  lattice-misfit  in¬ 
duced  strain  energy  is  added  to  the  Arrenius  law  to  accu¬ 
rately  account  for  the  effect  of  the  long-range  strain  field  on 
the  growth  patterns.  We  choose  the  substrate  to  be  GaAs,  and 
consequently,  the  lattice-misfit  strain  is  hydrostatic,  i.e., 

=  y* 8tj  with  y*  =  0.07.  The  nonzero  elastic  coefficients  for 
GaAs  (001)  (Ref.  26)  are  taken  to  be  Cn  =  118.8 
X  109  N/m2,  C12=53.8  X  109  N/m2,  and  C44=59.4 
X  109  N/m2,  when  the  growth  is  along  the  (001)  direction. 
For  the  substrate  with  growth  along  the  (111)  and  (113)  di¬ 


rections,  the  elastic  moduli  of  GaAs  (111)  and  GaAs  (113) 
are  transformed  from  GaAs  (001). 27  An  isotropic  substrate 
model,  called  Iso  (001),  is  also  introduced  for  the  purpose  of 
comparison  in  which  the  elastic  constants  of  GaAs  (001)  are 
replaced  with  the  corresponding  isotropic  constants.  In 
other  words,  in  the  isotropic  Iso  (001)  model,  the  elastic 
constants  C12  and  C44  are  the  same  as  those  of  GaAs  (001), 
i.e.,  C12=53.8  X  109  N/m2  and  C44=59.4  X  109  N/m2,  with 
Cn  =  172.6  X  109  N/m2  obtained  by  imposing  the  isotropic 
condition  C11-C12-2C44=0.3°  For  the  sake  of  easy  refer¬ 
ence,  all  the  elastic  coefficient  matrices  are  given  in  the  Ap¬ 
pendix. 

Once  we  have  defined  the  substrate,  we  then  use  our 
coupled  KMC  to  model  the  deposition  and  diffusion  pro¬ 
cesses.  Under  different  growth  parameters  (deposition  rate  F, 
temperature  F,  coverage  c,  and  interruption  time  tt)  we  simu¬ 
late  the  adatom  growth  behavior  over  a  grid  of  200  X  200  on 
the  surface  of  the  substrate.  We  remark  that  in  this  article, 
only  the  in-plane  or  the  so-called  two-dimensional  growth  is 
considered. 

We  further  point  out  that,  due  to  its  rapid  decay  behavior, 
the  strain  energy  field  is  not  evaluated  over  the  whole  prob¬ 
lem  domain  but  rather  only  over  a  circular  area  of  a  given 
radius  (e.g.,  thirty  grids25)  around  the  source  point.  This  de¬ 
cay  feature  can  be  observed  from  Fig.  1  where  the  elastic 
strain  energy  field  Fstr  on  the  surface  of  substrates  Iso  (001), 
GaAs  (001),  GaAs  (111),  and  GaAs  (113)  due  to  a  buried 
island  is  plotted.  Furthermore,  Fig.  1  also  shows  clearly  the 
direct  correlation  between  the  strain  energy  contour  shape 
and  the  substrate  orientation.  The  different  strain  energy  pat¬ 
terns  on  the  surface  of  Iso  (001),  GaAs  (001),  GaAs  (111), 
and  GaAs  (113)  will  result  in  different  island  distributions/ 
patterns,  as  we  will  discuss  below  using  our  coupled  KMC 
simulation. 

Shown  in  Figs.  2(a)-2(d)  are  the  island  orderings  and 
patterns  of  the  adatoms  on  the  surface  of  substrates  Iso  (001), 
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(a)  With  Estr  (Iso) 


(b)  WithEstr  (001) 


(d)  WithEstr  (113) 


(e)  Without  Estr 


(c)  With  E^r  (111) 


FIG.  2.  (Color  online)  Island  ordering  and  pattern  on 
the  surface  of  GaAs  with  strain  energy  EstT  of  Iso  (001) 

(a) ,  with  anisotropic  strain  energy  Estr  of  GaAs  (001) 

(b) ,  GaAs  (111)  (c),  and  GaAs  (113)  (d).  The  island 
ordering  and  pattern  on  the  surface  of  GaAs  without 
strain  energy  EstI  is  shown  in  (e)  for  comparison.  Dif¬ 
ferent  island  orderings  and  patterns  (red  dashed  lines) 
corresponding  to  the  GaAs  substrate  with  different 
growth  orientations  can  be  clearly  observed  from  (a)  to 
(d),  which  can  be  further  compared  to  the  elastic  strain 
energy  in  Fig.  1.  Other  simulation  parameters  are  T 
=750  K,  7=1.0  Ml/s,  c=20%,  and  ^=200  s,  on  a  grid 
of  200X200. 


GaAs  (001),  GaAs  (111),  and  GaAs  (113),  respectively.  The 
simulation  parameters  are  temperature  T=750  K,  flux  rate 
F=  1.0  Ml/s,  surface  coverage  c  =  20%,  and  interruption  time 
/)  =  200  s.  It  is  obvious  that  different  substrate  orientations 
result  in  different  growth  patterns,  due  mainly  to  the  different 
strain  energy  patterns  on  the  surface  of  the  substrate  (Fig.  1). 
Furthermore,  without  consideration  of  the  lattice-misfit  in¬ 
duced  strain  energy,  neither  uniform  size  nor  ordered  spatial 
ordering  can  be  observed  [Fig.  2(e)]  and  one  would  therefore 
expect  Ostwald  ripening.31,32 

In  order  to  show  the  island  distribution  more  clearly,  we 
plot  the  island  density  distribution  in  Fig.  3  where  the  num¬ 
ber  of  adatoms  per  unit  area  is  used  as  the  height.  From  Fig. 
3,  we  observe  again  that  without  strain  energy,  randomly 


(a)  With  Estr(lso)  (b)  With  E^fOOl)  (c)With  7^(111) 


(d)  With  Estdl  13)  (e)  Without  E^r 


FIG.  3.  (Color  online)  Three-dimensional  island  density  distribution  where 
the  height  is  proportional  to  the  number  of  adatoms  per  unit  area:  on  the 
surface  of  GaAs  with  strain  energy  7str  of  Iso  (001)  (a),  with  anisotropic 
strain  energy  7str  of  GaAs  (001)  (b),  GaAs  (111)  (c),  and  GaAs  (113)  (d). 
Different  island  orderings  and  patterns  (red  dashed  lines)  corresponding  to 
different  growth  orientations  can  be  clearly  observed,  which  can  be  further 
compared  to  the  elastic  strain  energy  in  Fig.  1.  Shown  in  (e)  is  the  island 
density  distribution  on  the  GaAs  substrate  without  strain  energy  7str.  Other 
simulation  parameters  are  7=750  K,  7=1.0  Ml/s,  c  =  20%,  and  ^=200  s,  on 
a  grid  of  200  X  200. 


ordered  islands  [Fig.  3(e)]  are  usually  formed  and  often  as¬ 
sociated  with  Ostwald  ripening;32  however,  with  the  strain 
energy  Estr,  isolated  and  ordered  islands  (for  both  size  and 
distribution)  can  be  observed  [Figs.  3(a)-3(d)].  Furthermore, 
different  island  patterns  are  correlated  closely  to  different 
growth  orientations  [GaAs  (001),  GaAs  (111),  and  GaAs 
(113)  to  Figs.  3(b)-3(d)],  with  the  latter  corresponding  to  the 
different  strain  energy  distributions  as  shown  in  Fig.  1.  These 
results  indicate  that  the  long-range  strain  energy  could  be 
one  of  the  very  important  factors  controlling  the  ordered  is¬ 
land  size  and  pattern.  Actually,  images  of  the  island  patterns 
on  different  anisotropic  substrates  from  experiments  have 
also  shown  that  a  direct  correlation  exists  between  the  island 
pattern  and  substrate  anisotropy  (or  growth  orientation),  as 
recently  presented  by  Jacobi33  where  InAs  QDs  were  grown 
on  GaAs  (001)  and  GaAs  (113)  substrates  and  by  Brune  et 
al.  ’  where  the  influence  of  the  elastic  strain  on  Ag  self¬ 
diffusion  and  nucleation  was  investigated. 

III.  EFFECT  OF  DIFFERENT  GROWTH  PARAMETERS 

In  our  previous  paper,25  we  studied  the  effect  of  growth 
parameters  (temperature,  flux  rate,  coverage,  and  interruption 
time)  and  showed  some  typical  island  patterns  for  different 
growth  parameters.  In  that  paper,  however,  no  discussion  is 
given  on  the  effect  of  substrate  anisotropy  (or  orientation)  on 
the  growth  pattern.  Therefore,  in  what  follows,  we  investi¬ 
gate  the  effect  of  the  growth  parameters  on  the  number  of 
islands  and  average  size  of  the  islands.  Again,  we  use  the 
four  different  substrates  in  our  simulation:  the  isotropic  sub¬ 
strate  Iso  (001)  and  three  anisotropic  GaAs  substrates  with 
different  growth  directions,  namely,  GaAs  (001),  GaAs 
(111),  and  GaAs  (113).  For  all  the  examples  below,  the  grid 
size  is  200  X  200  and  the  strain  energy  is  included.  We  dis¬ 
cuss  our  results  for  the  temperature,  coverage,  and  interrup¬ 
tion  time  separately.  The  effect  of  flux  rate  on  the  island  size 
was  discussed  in  our  previous  paper  where  it  was  identified 
that  a  small  size  island  distribution  usually  corresponds  to  a 
high  flux  rate  and  a  large  size  island  distribution  to  a  low  flux 
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FIG.  4.  (Color  online)  Number  of  islands  and  average  island  size  vs  tem¬ 
perature  for  island  growth  on  substrates  with  different  orientations  for  short 
temperature  range  (550-750  K).  Other  simulation  parameters  are  fixed:  7 
=  1.0  Ml/s,  c=20%,  and  fi=200  s. 
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FIG.  6.  (Color  online)  Number  of  islands  and  average  island  size  vs  cover¬ 
age  for  island  growth  on  substrates  with  different  orientations.  Other  simu¬ 
lation  parameters  are  fixed:  7=750  K,  7=1.0  Ml/s,  and  ^=200  s. 


Figures  4  and  5  show  the  variation  of  the  number  of 
islands  and  average  island  size  versus  temperature.  Other 
fixed  growth  parameters  are  flux  rate  F=l  Ml/s,  coverage 
c  =  20%,  and  interruption  time  f;  =  200  s.  We  observed  from 
Fig.  4  that,  under  given  growth  conditions,  the  number  of 
islands  and  average  island  size  are  insensitive  to  the  substrate 
anisotropy  (or  growth  orientation).  For  all  the  substrates, 
with  increasing  temperature  to  750  K,  the  island  size  in¬ 
creases  and  the  number  of  islands  decrease.  In  our  previous 
discussion,  we  further  identified  by  looking  at  the  island  size 
ordering  and  distribution  that  the  optimal  temperature  was 
around  750  K.25  The  island  distribution  could  be  still  irregu¬ 
lar  if  temperature  is  too  low  (below  this  critical  value)  and  be 
unstable  if  temperature  is  too  high  (above  this  critical 
value).25  Therefore,  only  around  the  optimized  temperature 
(i.e.,  T=1 50  K)  can  one  clearly  observe  ordered  and  regular 
island  distribution.  To  further  confirm  our  previous  observa¬ 
tion,  we  show  here  also  the  variation  of  the  average  island 
size  and  number  of  islands  versus  temperature  for  a  large 
temperature  range  [Fig.  5  for  GaAs  (001)  only].  It  is  ob- 


Temperature  (K) 

FIG.  5.  (Color  online)  Number  of  islands  and  average  island  size  vs  tem¬ 
perature  for  island  growth  on  substrate  GaAs  (001)  for  long  temperature 
range  (550-950  K).  Other  simulation  parameters  are  fixed:  F=  1.0  Ml/s, 
c=20%,  and  tt= 200  s. 


served  from  Fig.  5  that  the  average  size  of  the  island  reaches 
a  maximum  at  T=1 50  K  where  the  number  of  islands 
reaches  a  minimum. 

Figure  6  shows  the  variation  of  the  number  of  islands 
and  average  island  size  on  the  surface  of  substrates  Iso  (001), 
GaAs  (001),  GaAs  (111),  and  GaAs  (113)  for  different  cov¬ 
erages  c.  The  other  fixed  parameters  are  temperature  T 
=  750  K,  flux  rate  F=l  Ml/s,  and  interruption  time  tt 
=  200  s.  It  is  interesting  that  with  increasing  coverage,  while 
the  number  of  islands  in  the  grid  is  nearly  a  constant  value 
(around  50),  the  average  island  size  increases,  with  GaAs 
(113)  having  the  largest  value  followed  by  GaAs  (111).  The 
average  island  size  on  the  surface  of  the  substrate  GaAs 
(001)  has  the  smallest  value. 

Finally,  Fig.  7  shows  the  effect  of  the  growth  interrup¬ 
tion  time  tt  on  the  number  of  islands  and  average  island  size. 
The  fixed  growth  parameters  are  T=1 50  K,  F=  1.0  Ml/s,  and 
c  =  20%.  It  is  observed  clearly  from  Fig.  7  that,  with  increas¬ 
ing  interruption  time,  the  number  of  islands  decreases  while 
the  average  island  size  increases.  Furthermore,  at  large  inter¬ 
ruption  time  (f;=200  s  for  our  case),  while  the  number  of 
islands  is  nearly  the  same  for  different  substrate  orientations 
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FIG.  7.  Number  of  islands  and  average  island  size  vs  interruption  time  for 
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(or  different  material  anisotropies),  different  substrate  orien¬ 
tations  predict  slightly  different  island  sizes  with  GaAs  (111) 
having  the  smallest  size  among  ah  the  substrates. 

IV.  CONCLUSIONS 

In  this  paper,  a  coupled  kinetic  Monte  Carlo  algorithm  is 
proposed  to  simulate  the  growth  of  adatom  islands  on  the 
surface  of  isotropic  and  anisotropic  substrates.  The  substrates 
considered  are  isotropic  Iso  (001)  and  anisotropic  substrates 
GaAs  (001),  GaAs  (111),  and  GaAs  (113).  The  unique  fea¬ 
ture  associated  with  our  coupled  KMC  is  the  correct  incor¬ 
poration  of  the  lattice-misfit  induced  strain  energy  into  the 
diffusion  process  governed  by  the  Arrenius  law.  Our  simula¬ 
tion  results  are  consistent  with  previous  experimental  obser¬ 
vations  which  show  clearly  that  the  growth  parameters  (tem¬ 
perature,  coverage,  and  interruption  time)  are  important 
parameters  in  controlling  the  island  size/distribution.  Our  re¬ 
sults  also  show  that  without  considering  the  strain  energy, 
only  randomly  ordered  islands  can  be  formed,  which  is  con¬ 
trolled  by  the  nearest  and  next  nearest  binding  energy  among 
the  adatoms  and  substrate  atoms.  When  the  strain  energy 
held  is  included,  however,  isolated  and  ordered  island  pat¬ 
terns  (both  in  size  and  distribution)  can  be  observed.  Further¬ 
more,  our  simulation  for  self-organized  island  growth  on  the 
surface  of  the  isotropic  Iso  (001)  and  anisotropic  GaAs 
(001),  GaAs  (111),  and  GaAs  (113)  indicates  that  the  growth 
pattern  is  closely  correlated  to  the  substrate  anisotropy  (or 
substrate  orientation). 
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APPENDIX:  ELASTIC  COEFFICIENTS 
OF  SUBSTRATES  ISO  (001),  GAAS  (001), 

GAAS  (111),  AND  GAAS  (113) 

An  elastic  isotropic  model,  called  Iso  (001)  is  used  for 
the  purpose  of  comparison  in  which  the  elastic  coefficients  of 
GaAs  (001)  are  made  isotropic  (C11-C12-2C44=0).30  For 
the  GaAs  (111)  model,  the  substrate  coordinates  are  such  that 
the  x  axis  is  along  [11-2],  y  axis  along  [-110],  and  z  axis 
along  [111],  while  for  the  GaAs  (113)  model,  the  v  axis  is 
along  [33-2],  y  axis  along  [-110],  and  z  axis  along  [113]. 
The  global  material  properties  of  GaAs  (111)  and  GaAs 
(113)  are  obtained  by  the  well-known  tensor  transformation34 
from  GaAs  (001).  The  inclined  material  property  is  obtained 
by  two  transformations.  For  GaAs  (111),  we  first  anticlock¬ 
wise  rotate  xy  plane  around  z  axis  by  7t/4  (view  from  posi¬ 
tive  z  axis);  then  we  anticlockwise  rotate  the  new  xz-plane 
around  the  new  y  axis  by  7t/2-cos_1(V2/3)  (view  from  posi¬ 
tive  new  y  axis).  Similarly,  for  GaAs  (113),  we  first  anti¬ 
clockwise  rotate  xy  plane  around  z  axis  by  7t/4  (view  from 
positive  z  axis);  then  we  anticlockwise _rotate  the  new  yz 
plane  around  the  new  y  axis  by  tan-1  (3/  ^2)  (view  from  posi¬ 
tive  new  y  axis).  The  elastic  coefficients  thus  obtained,  along 
with  those  for  the  Iso  (001)  and  GaAs  (001),  are  listed  be¬ 
low:  nonzero  elastic  coefficients  of  Iso  (001)  (X109  N/m2), 


Cn  =  172.6;  C12  =  53.8;  C44  =  59.4;  (Al) 

nonzero  elastic  coefficients  of  GaAs  (001)  (X109  N/m2), 
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